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ABSTRACT 


Objective: The objective is to obtain static and quasi-static deformation of a uniform half-space due to a center of rotation. 

Methodology and Results: The Galerkin vector approach has been used to calculate deformation field at an arbitrary point 
of an elastic half-space. Closed form analytical expressions for the displacements and stresses caused by a center of rotation 
buried in a homogenous, isotropic, perfectly elastic half-space are derived. The quasi-static deformation field for a viscoelastic 
medium has been obtained by applying the correspondence principle of linear viscoelasticity to the associated elastic solution. 
Explicit expressions giving the quasi-static deformation of a uniform half-space caused by a center of rotation are obtained when 
the medium is elastic in dilatation and Kelvin, Maxwell or SLS type viscoelastic in distortion. 

Conclusion: The explicit expressions for the displacements and stresses in an elastic and viscoelastic medium due to a center 
of rotation source have been obtained. Numerical results are shown graphically for displacements and stresses. 
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INTRODUCTION 

Nuclei of strain are the concentrated sources of a displace¬ 
ment field and are built up from the simple superposition of 
single forces which are acting at a point in the elastic medium 
(Love 1 ). Therefore, the displacement field owing to nuclei of 
strain is necessary for applications to crustal deformation. 
Analytical expressions for three-dimensional static displace¬ 
ment fields owing to a step-type single force in an elastic 
half-space, have been obtained by Mindlin 2 by using Galer- 
kin’s method. Mindlin and Cheng 3 derived the solutions in 
the form of the Galerkin vector for various nuclei of strain 
in a uniform half-space by the process of superposition, dif¬ 
ferentiation and integration. Many theoretical formulations 
have been developed (Okada 4 ), which were describing the 
deformation of an isotropic, homogeneous, and semi-infinite 
medium. Analytical expressions for the quasi-static surface 
displacements due to a vertical strike-slip fault in a Kelvin or 
Maxwell viscoelastic half-space were determined by Singh 
and Rosenman 5 by applying the correspondence principle 
of linear viscoelasticity. The correspondence principle has 
been broadly used to calculate the quasi-static deformation 


of a viscoelastic half-space by a point or extended sources 
(see. e.g. [6-10]). Singh and Singh 11 have identified the com¬ 
binations in which moduli occur in the expressions for the 
displacements, strains and stresses in a uniform elastic half¬ 
space due to buried sources. Using analytical integration, the 
displacement field in two welded elastic half-spaces due to a 
finite rectangular fault has been obtained by Singh et al 12 and 
they have also compared it with the corresponding field in an 
elastic half-space and in an infinite medium. To model the 
ground defonnation in volcanic areas, Singh et al 13 used four 
axially-symmetric source models in an elastic half-space and 
also compared it with the corresponding field due to a center 
of dilatation. 

Cochard et al 14 showed that the observations of seismic rota¬ 
tional motions will give important new information referring 
to the Earth’s surface and are complementary to those found 
from the observations of the translational motions of the 
Earth’s surface using conventional seismometers. Cowsik et 
al 15 resolved that to detect rotational motions, the basic de¬ 
sign concept of using a torsion balance as a filter is validated 
and can be implemented for the construction of rotational 
seismometers. 
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The static displacements from the seismic recordings and 
identifying translation signals caused by rotation can be 
estimated by using rotational motions (Trifunac and Todor- 
ovska 16 ). Rotational seismology is of great interest to a wide 
range of geophysical disciplines, including strong-motion 
seismology, seismic tectonics, earthquake engineering, and 
geodesy as well as to physicists using Earth-based observa¬ 
tories for detecting gravitational waves generated by astro¬ 
nomical sources (Lee et al 17 ). The solutions for the displace¬ 
ment field produced by a center of rotation is also useful for 
many purposes and in bio-mechanical research. 


RESEARCH METHODLOGY 

In this paper, we study the 3-D deformation of a uniform 
half-space caused by a center of rotation by using the Galer- 
kin vector approach. Explicit expressions for the static 
strains can be easily obtained with the help of strain -dis¬ 
placement relations and the stresses follow immediately by 
using stress-strain relations. The correspondence principle of 
linear viscoelasticity has been used to obtain the quasi-static 
displacements, strains and stresses. 

The paper has been divided into two parts. Part-A deals with 
the static deformation field while the quasi-static deforma¬ 
tion field is considered in Part-B. 


PART- A: STATIC DEFORMATION FIELD 
THE GALERKIN VECTOR 

Mindlin and Cheng 3 expressed the displacement vector u in 
term of the Galerkin vector (} through the relation 

2pu = 2(l - a)v 2 G- v(v.g) (2.1) 

where p is the shear modulus and o is the Poisson’s ratio. 
From equation (2.1), the displacement u [ in the x. - direction 
is given by 


2(l - o)V 2 G; --^(v.g) 

CX • 


( 2 . 2 ) 


where G = G i e i and e i denotes the unit vector in the x. - 
direction. 

The strain-displacement relations are 


e 


ij 


1 

2 


3 u i 3uj 

dx ■ dx- 

< 1 > 


where e , f is the strain tensor. 


(2.3) 


The stress-strain relations are 

*ij=M i j u n,n +1 Ve ij ( 2 . 4 ) 

where r //' is the stress-tensor, 8 „ is the Kronecker delta and 
4,// are Lame’s constants. From equation (2.1) and (2.3), 

= K, n = divu = (v.g) (2.5) 

2 // 

The expressions for the Galerkin vectors for a vertical and 
horizontal concentrated force in a homogeneous, isotropic, 
elastic half-space were given by Mindlin and Cheng >. Let 
the uniform half-space be z > 0 with a stress-free boundary 
at z=0 and the point source is located at the point (0,0, c) of 
the upper half space as shown in Fig. 1. R, is the distance 
between the observer at (x, y, z) and the source at (0,0, c); 
R 2 is the distance between the observer at (x, y, z) and the 
image at (0,0, - c), 

where 

R~ = x 2 + y 1 +(z-cf 
R 2 2 = x 2 + y 2 +(z + cf ( 2 . 6 ) 

CENTER OF ROTATION 

Combining two coplanar single couples with moments equal 
in magnitude and sense, we get a torque or center of rotation, 
in which the net force is zero, whereas the net moment is 
non-zero, as shown in Fig.2. 

For further study, we use (x, y, z) co-ordinate instead of 
(x l ,x 2 ,F 0 ). 

Taking strength of center of rotation along x-axis, y-axis and 
z-axis as F„ which acts at the point (0,0,c), the expressions for 
Galerkin vectors are given as under: 

Center of rotation about x-axis 


G = 


8/r(l - cr) 


z — c z —3c 2 c 2 (z + c) 

- +-+-- r —- + 


R\ R2 


R,- 


4 (l - cr) (l - 2cr) log (R 2 + z + c) 
1 + 4cr(l-2o-) + 3 2c 2 


K 




R, 


R 2 


R 2 ' 


4(1- 2 < 7 ) [(1 - cr) z - ccr] 


R 2 (R 2 +z + c) 


(2.7) 
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Center of rotation about y-axis 

z-c z-3e 2 c 2 (z + c) 

R, R 2 R 2 1 
4(l-(r)(l-2c7)log(i? 2 + z + c) 
1 4<t(1-2ct) + 3 2 c 2 


G = 


8;r(l - cr) 


4(l-2<r)[(l-tr)z-trc] 


R 2 ( R 2 + z + c 

Center of rotation about z-axis 


Center of rotation about z-axis 

Using equations (2.2) and (2.9), we obtain the following ex¬ 
pressions for the displacement components for a center of 
rotation in the z- direction 


,( 3 )_ F °y 

1 Ann 


( 2 . 8 ) 


(3)_ 


Fx 


Ann 


J_ 1 

r^ + rJ 


J_1_ 


G = 


F 


8^(1-cr) 


g i>') 


1 1 2c 

-1-1- — 

R, id 7? 3 


4(1 - a) (1-2a) 


7? 2 + z + c 


u® = 0 


(3.7) 

(3.8) 

(3.9) 


(2.9) 


DISPLACEMENT FIELD 

Center of rotation about x-axis 

Using equations (2.2) and (2.7), we obtain the following ex¬ 
pressions for the displacement components for a center of 
rotation in the x - direction 


, w = F o*y 

Ann 


2(l-2cr)(27? ; , + z-i-c) 6z 
R 2 3 (R 2 +z + c) 2 R 2 5 


STRAIN AND STRESS FIELD 

The strain components e can be easily obtained by using the 
displacements components u i as expressed in previous sec¬ 
tion in relation (2.3) and the stresses follow from equation 
(2.4), which are given as below: 

STRESSES 

Stress components due to center of rotation about x-axis 


6c 6(1-2 o)(z+c) 2(l-2<r)(2R, + z + c) 

R 2 r i + R 2 (R 2 + z + c) 2 


(3.1) 


F_ 

Ann 


in _ F oy 


z-c z — c 2(1-2ct) 

R 3 R 2 3 R 2 (F, + z 4- c) 


+r 


R' 


Ann 


2(l-2cr)(2F, +z + c) 6 z 
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1 6z(z + c) 3-4cr 

'W ~r 5 


, w _ F °y 

T " ~2k 
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(3.2) 


(3.3) 


4(l-2o)(2R,+z + c) 2 

2(1-2ct)(z + c) 

R 2 5 + Z + c) 3 

R 2 (R 2 +z + c ) 2 

30 z 


+ R 2 





(4.1) 


Center of rotation about y-axis 


Using equations (2.2) and (2.8), we obtain the following ex¬ 
pressions for the displacement components for a center of 
rotation in the y- direction 


T w-Tr 

22 “ 2/r 


,(>) = . 


Ann 


z — c z—c 2(1 — 2<r) 

R j" R 2 R 2 (F, + z + c) 


+JC‘ 


2(\-2<j)(2R 2 +z + c) 6 z 
R 2 3 (R 2 + z + c ) 2 R 2 


(2) _ F 0 xy 


u^ ] = 


Ann 

Fx 


Ann 


2(l — 2 <t)( 27?-, + z + c) 6z 
R 2 (R 2 +z + cf R 2 

1 6z(z + c) 3-4cr 


R; 


R 
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(3.4) 

(3.5) 

(3.6) 


3(z-c) 3(5 z-c-4ct(z + c)] 

‘Ur/ r? 

6(l-2tj)(2R 2 +z + c) 


r/(r 2 


4(l-2ff)(2R,+ z + c) 2 

2(l-2a)(z+c) 

R/(R,+z + c) 3 

R, s (R,+z + e) 2 
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3 ( 

z-c) 3(3 z-c) _ 4(1 - 2o)(lR 2 + z + c) 


R l R 2 r 2 3 ( r 2 + z + cf 



60z 4(l-2a)(z + c) 


F„x 


R 2 R 2 5 ( R 2 +Z + cf 


An 
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8(l-2a)(2« 2 +z + c) 2 



R 2 ( R 2 + z + c) 3 







, m _ F o y 
12 ~ 4k 


(4.4) 


3(z-c) 3(3z-c) i 4(l-2a)(2it,+z + c) 
R i + fC(ft, + z + c) 2 

60x 3 4(l-2a)(z + c) 

^“^(R+z + c) 2 
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( 2 ) = 3F„ 
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20z(z + c) 1 1 
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*, 5 
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+x 


20z(z + c) i 1 1 
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(4.12) 

(4-6) Stress components due to center of rotation about z-axis 


Stress components due to center of rotation 
about y-axis 


3 [z-c) 3(5z-c-4(t(z + c)) 

R* 

6(1 — 2<j)(2i? 2 + z + c) 
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(4.14) 

(4.15) 

(4.16) 

(4.17) 

(4.18) 


PART B. QUASI-STATIC DEFORMATION 

BASIC EQUATIONS 

The stress-strain relations (2.4) for an isotropic, linear elastic 
medium, may be expressed in the form 


Pj = 2/i£ iJ -e=3Ki} 


(5.1) 
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where 

0=T n + T 22 + T 33 

4 — e n + e 22 + e 33 


Pij -deviatoric stress =*q-- 08 ij , 

£„ = deviatoric strain =e , and 

y y 2 lJ 


2 

v= A + -« is known as bulk modulus 
3 

The stress-strain relations for an isotropic, linear viscoelastic 
medium, may be expressed in the form 


F o 

4k 


„ T z-c z-c 6zy~ 
2J,(t)\—^ + 


«( 

-64(0 


R 3 R 2 3 R 2 

1 y 2 (2R 2 +z + c) 


R 2 (R 2 + z + c) R 3 [R^ + z + c ) 2 


(5.4) 


, (i) _ - F *y 


4 n 


2j\t) 

, 64(0 

F? 


1 6 z(z + c) 1 
~^ + R. 5 + F : 


vR 


(5.5) 


Quasi-static displacements due to center of rotation 
about y-axis 


Pu ( s ) = 2 P{ s ) £ j i s ). 0( s ) = 3r ( s ) &i s ) ( 5 - 2 ) 


where 


Mi*) 


QM 

2 P,(s) 


k(s) 


QM 

3 P 2 (s) 


where are the polynomials of degree (positive in¬ 

tegers), respectively. 

From equations (5.1) and (5.2), we note that on replacing the 
elastic moduli P and K by the transformed modidi /7 and 
^respectively, the Laplace transformed viscoelastic solution 
can be obtained from the corresponding transformed elastic 
solution. This is the well-known correspondence principle of 
viscoelasticity. 


F o 

4k 


- , n. i z — c z — c 6 zx 2 


-«,(<) 


1 x 2 (2 R 2 + z + c) 
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6 Q 3 (t)(2R 2 + z + c) _ 12 J x (t)i 
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-Fx 
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- , 1 6z(z + c) 1 
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(5.8) 


Quasi-static displacements due to center of rotation 
about z-axis 


The elastic moduli // and X occur in various combinations 
in the expressions for the displacements, strains, and stresses. 
For these combinations, correspondence between the static 
solution for an elastic medium and quasi-static solution in the 
time domain for a viscoelastic medium was given in table- 
1 (Singh and Singh 10 ). We can easily obtain the quasi-static re¬ 
sponse of a viscoelastic half-space to various sources from the 
static responses given in Part-1 by using table-1. 

The quasi-static displacements, strains and stresses for center 
of rotation are obtained when the medium is elastic in dilata¬ 
tion and kelvin. Maxwell or SLS type viscoelastic in distor¬ 
tion. 

QUASI-STATIC DISPLACEMENTS 

Quasi-static displacements due to center of rotation about x- 
axis 
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6Q i (t)(2R 2 +z + c) 
R 2 (. R 2 +z + cy 
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QUASI-STATIC STRESSES 

Quasi-static stresses due to center of rotation about x-axis 
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Quasi-static stresses due to center of rotation about y-axis 
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f 20 z 
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K R2 


w = 


D(l-e(-T)) 

3 

(d 2 + 


(6.4) 


where t = — and ^ is relaxation time for Kelvin model. 


Maxwell model 


(5.23) 


Quasi-static stresses due to center of rotation about z-axis 
are same as static expressions of stresses given in section 
4 . 1 . 3 . 

NUMERICAL RESULTS 

We define dimensionless epicentral distance (D), dimension¬ 
less horizontal displacement (U), dimensionless vertical dis¬ 
placement (W), dimensionless stresses ( T x , ) by the rela¬ 

tions 

d= 1,u = ^ w= A u ®' Tx = a tJ 1) and t, = A T n (3) where 

a-- i‘ v " is dimensionless constant. 

F « 

For further calculations, the following assumptions are con¬ 
sidered: the half-space z > 0 is stress-free at z=0, the bulk 
modulus K and the shear modulus A of viscoelastic half- 

/£* 5 

space satisfy Poisson’s condition - = -. We also assume 
that t>0 and the source time function is the Heaviside step 
function. Using the values of various auxiliary functions are 
given by Singh and Singhll, considering the material as an 
elastic in dilatation and Kelvin, Maxwell, or SLS viscoelas¬ 
tic in distortion. 
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1 +T 3 [ 6 
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z)(i+r) 
(z? 2 + 1 )^ 


(u 2 +i) 5 ^ z) 2 +1 j 2+1 j 


(6.5) 
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where r = — and t-, is relaxation time for this model. 


SLS model 
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( 6 . 8 ) 


STATIC DISPLACEMENTS 

-D 

3 lix-ici^ip 1 + 1 ) 5+1 j 

From equations (3.1) -(3.9), the dimensionless static displace¬ 
ments are as follows: 

(A+\j 

H ((z> 2 4 + ij 


u=- 


1 

1 : (1 - 2cr) 

(z> 2 + 1 ) 2 

( Z) “ + 1 ) j^(Z) 2 +l) 2 +lj 


where T -~r and U is relaxation time for SLS model. 

‘1 z 

STATIC STRESSES 

From equations (4.1) -(4.18), the dimensionless static stress¬ 
es are as follows: 


(6.9) 

( 6 . 10 ) 


r=- 


-3 D 1 


( 6 . 1 ) 


W = - 


D 


[p 2 + 1) 2 

QUASI-STATIC STRESSES 

From equations (5.12) -(5.17), the dimensionless quasi-static 
stresses for various viscoelastic models are as follows: 


(Zf+l) 


(6.2) Kelvin model 


QUASI-STATIC DISPLACEMENTS 

From equations (5.3) -(5.11), the dimensionless Quasi-static 
displacements for various viscoelastic materials is obtained 
as under: 

Kelvin model 
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SLS model 

-M [ I 


DISCUSSION 

Fig. 3. shows the effect of the Poisson’s ratio on the variation 
of the dimensionless horizontal displacement (U) with the di¬ 
mensionless epicentral distance (D). We note that, horizontal 
displacement increases with the decrease of Poisson’s ratio 
and epicentral distance. It vanishes at D=1.272 for <7 = 1. 

Fig. 4. shows the variation of the dimensionless vertical dis¬ 
placement (W) with the dimensionless epicentral distance. We 
note that, it increases with epicentral distance and assumes 
maximum value W=0.35 at D=1 after that it decreases rapidly 
as epicentral distance increases. 

Fig. 5. shows the variation of the dimensionless stress L with 
the dimensionless epicentral distance for various values of <T 
. We note that, it increases with the increase of Poisson’s ratio 
and assumes minimum value -0.8 at =1/2 for D= 0.6 after 
that it increases with epicentral distance and becomes constant 
as D tends to infinity. 

Fig. 6. shows the variation of the dimensionless stress T z with 
the dimensionless epicentral distance. We note that, stress T, 
decreases with epicentral distance and assumes minimum 
value at D=1 after that it increases with epicentral distance. 

Fig. 7. shows the variation of horizontal displacement for the 
three models namely Kelvin, Maxwell and SLS with the epi¬ 
central distance for different values of T. We note that in case 
of Kelvin model, horizontal displacement decreases with the 
epicentral distance as shown in Fig.7(a) but in case of Max¬ 
well model(Fig.7(b)) and SLS model(Fig.7(c)), it decreases 
first and then increases with epicentral distance. 

Fig. 8 shows the variation of the vertical displacement for the 
three models namely Kelvin, Maxwell and SLS with the epi¬ 
central distance for various values of T as shown in Fig. 8(a), 
Fig. 8(b) and Fig. 8(c) respectively. The maximum and mini¬ 
mum values of the vertical displacement are given in Table 1. 
Vertical displacement assumes the maximum value in case of 
Maxwell model and minimum in case of Kelvin model for a 
particular value of time. 

Fig. 9. shows the variation of the dimensionless stress T x for 
the three models namely Kelvin, Maxwell and SLS with the 
epicentral distance for various values of T as shown in Fig. 
9(a), Fig. 9(b), and Fig. 9(c), respectively. We note that in case 
of Kelvin(Fig.9(a)) stress is same for various values of time 
except at T=0. In all the three cases, stress decreases first and 
then increases with the epicentral distance. 


CONCLUSIONS 

The explicit expressions for the displacements and stresses 
in an elastic and viscoelastic half-space due to a center of 
rotation source have been obtained by using Galerkin vec¬ 
tor approach. Numerical results are shown graphically for 
displacements and stresses for different values of Poisson’s 
ratio and three viscoelastic models namely Kelvin, Maxwell 
and SLS. 
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Table i: Below table containsthe maximum and mini¬ 
mum values of the vertical displacement. 


Time(T) 

Max.(Min.) value of vertical displacement (W) 


Kelvin 

Maxwell 

SLS 

o 

o(o) 

0.3(0) 

0.35(0) 

1 

0.2(-0.2l) 

0.65(0) 

0.49(0) 

2 

o. 3 (-o. 3 ) 

1(0) 

0.55(0) 

3 

o- 3 i(-c>- 3 i) 

1.4(0) 

0.62(0) 



Figure 3: Variation of the dimensionless horizontal displace¬ 
ment (U) with the dimensionless epicentral distance (D) for 
various values of (7 . 



Figure 1 : Geometry of a point source in a uniform half-space. 


( 32 )-( 23 ) 






EPICENTRAL DISTANCE(D) 


Figure 5: Variation of the dimensionless stress t x with the di¬ 
mensionless epicentral distance for various values of a. 


Figure 2: Three fundamental torques, r 1 -T 3 ,r 3 -T 3 ,T 2 -T 1 
are center of rotations about the and X 3axes respectively. 
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EPICENTRAL DISTANCE 


Figure 6: Variation of the dimensionless stress T z with the 
dimensionless epicentral distance. 



EPICENTRAL DISTANCE(D) 


Figure 7(a): Variation of the dimensionless horizontal dis¬ 
placement with the dimensionless epicentral distance for the 
various values of T for kelvin model. 



EPICENTRAL DISTANCE(D) 


Figure 8(a): Variation of the dimensionless vertical displace¬ 
ment with the dimensionless epicentral distance for various 
values of T for kelvin model. 



Figure 8(b): Variation of the dimensionless vertical displace¬ 
ment with the dimensionless epicentral distance for various 
values of T for Maxwell model. 



EPICENTRAL DISTANCE(D) 


Figure 7(b): Variation of the dimensionless horizontal dis¬ 
placement with the dimensionless epicentral distance for vari¬ 
ous values of T for Maxwell model. 



EPICENTRAL DISTANCE(D) 


Figure 7(c): Variation of the dimensionless horizontal dis¬ 
placement with the dimensionless epicentral distance for vari¬ 
ous values of T for SLS model. 



Figure 8(c): Variation of the dimensionless vertical displace¬ 
ment with the dimensionless epicentral distance for various 
values of T for SLS model. 



EPICENTRAL DISTANCE(D) 


Figure 9(a): Variation of the dimensionless stress T x with the 
dimensionless epicentral distance for various values of T for 
kelvin model. 
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EPICENTRAL DISTANCE(D) 


Figure 9(b): Variation of the dimensionless stress ^ with the 
dimensionless epicentral distance for various values of T for 
Maxwell model. 



Figure 9(c): Variation of the dimensionless stress T x with the 
dimensionless epicentral distance for various values of T for 
SLS model. 
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